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ABSTRACT 

Type la supernovae (SNe la) occur in both old, passive galaxies and active, star-forming galaxies. 
This fact, coupled with the strong dependence of SN la rate on star formation rate, suggests that SNe 
la form from stars with a wide range of ages. Here we show that the rate of SN la explosions is about 
1% of the stellar death rate, independent of star formation history. The dependence of SN la rate on 
star formation rate implies a delay time distribution proportional to t-0-5±0-2^ j^^ie single degenerate 
channel for SNe la can be made to match the observed SN la rate - SFR relation, but only if white 
dwarfs are converted to SNe la with uniform efficiency of '^1%, independent of mass. Since low-mass 
progenitors are expected to have lower conversion efhciencies than high mass progenitors, we conclude 
that some other progenitor scenario must be invoked to explain some, or perhaps all, SNe la. 
Subject headings: supernovae: general — stars: evolution 



1. INTRODUCTION 

Type la supernovae (SNe la) are, next to gamma ray 
bursts, among the most powerful transient objects in 
the Universe. SNe la are widely believed to be ex- 
plosions of C-t-0 white dwarfs (WDs) whose mass has 
grow n by accretion to the Chandrasekhar mass of 1.4 
M0 (jChandrasekhail Il931h . Evidence for this includes 
total energy released (consistent with conversion of C-l-0 
to Fe), lack of hydrogen in the maximum light spectra 
(consistent with an object whose envelope has been lost), 
light curve shape (broadly in agreement with energy de- 
position from radioactive decay of Fe-peak elements, as 
expected for exploding C-l-0 WDs), and the occurrence 
of SNe la in old (~10 Gyr) stellar populations containing 
only low mass stars (for which exploding WDs represent 
the natural, and perhaps only possible, SN la progeni- 
tor). 

Beyond this consensus on the exploding WD model, 
the precise nature of the progenitors of SNe la remains 
poorly constrained. Progenitor models are broadly clas- 
sified as "single-degenerate" (SD), in which a white 
dwarf grows in m ass by accretion from an evolving bi- 
nary companion (|Whelan fc IbenI [TOTSh . and "double 
degenerate" (DP), in which two wh ite dwarfs merge 
(|Iben fc Tutukovlll984l : IWebbinkl[l98l . A variety of evo- 
lutionary pat hs may lead to SD and DD events (e.g . 
Livid "2001: Yu ngelsonI [2001 iParthasarathv et all l2007l : 
GreggiQ.2005.) . 

The fact that SNe la are observed in elliptical galax- 
ies has led to the persistent and widespread presump- 
tion that SN la progenitors are old, low-mass stars. 
Yet over the years there have been several hints sug- 
gesting that t he SN la rate is enhanced in star-forming 
galaxies (e.g. lOemler fc TinslevI Il979t Ivan den Berghl 

^ Department of Physics and Astronomy, University of Victo- 
ria, PO Box 3055, Stn CSC, Victoria, BC V8W 3P6, Canada; 
pritchet@uvic.ca 

^ Department of Astronomy, University of Toronto, 50 
St. George St., Toronto, ON M5S 3H4, Canada; how- 
ell@astro.utoronto.ca 

^ Department of Physics (Astrophysics), University of Oxford, 
Denys Wilkinson Building, Keble Road, Oxford, OXl 3RH, United 
Kingdom; suIIivan@astro.ox.ac.uk 



119901: iDella Valle fc Li\dol I1994D. Two re cent papers 
(jMannucci et all l2005t ISullivan et all l2006f ) have mar- 
shalled evidence for a strong dependence of the SN la 
rate on star formation rate (SFR), thu s linking SNe la 
with y oung, massive stellar populations. iMannucci et al.l 
((2005h made use of a sample of low reds hift super- 
novae discovered in 5 targeted searches (ICappellaro et all 
Il999( l. and demonstrated a strong dependence of spe- 
cific SN la rate (rate per unit mass) on morphology 
or galaxy colour. This r esult was confirmed and ex- 
tended by ISullivan et al.l (|2006( ). who analyzed a sam- 
ple of spectroscopically-confirmed SNe la (redshift z=0.2- 
0.75, spanning cosmic time 7-11 Gyr after the Big Bang) 
from the Supernova Legacy Survey (jAstier et al.ll2006l) 
(SNLS). 

The SN la rat e - SFR correlation found by 
ISullivan et al.l (|2006f ) is shown in the lower panels of Fig. 
1 (both axes are normalized by galaxy mass). Passive 
galaxies are plotted as a single data point in the lower 
left hand panel; most of these galaxies are old. The vari- 
ation of SN la rate from passive galaxies to the most 
active galaxies is dramatic: passive galaxies possess a 
mass-normalized SN la rate that is a factor of ten or 
more lower than found for active star-bursting systems. 

In this paper we present a simple model for the SN 
la rate - SFR relation in terms of stellar evolutionary 
timescales. This model, which is based on the SD sce- 
nario, allows an estimation of the efficiency of conversion 
of white dwarfs into SNe la. This provides, as we shall 
see, a constraint on the nature of SN la progenitors. 



2. THE MODEL 



2.1. 



White Dwarf Formation Rates 

A first step towards understanding the rate of SNe la 
in galaxies is a calculation of the rate at which WDs 
are formed in different stellar populations. It is well- 
known that WDs (M < 1.4 Mq) form as the end-point 
of stellar evolution for stars wit h initial masse s M < 8 
Mq (e.g. [Wcidemann 2000; Kah rai et ani2007l . and ref- 
erences therein). The rate at which white dwarfs form 
clearly depends on evolutionary time scales as a function 
of initial stellar mass, and the initial mass function or 
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Fig. 1. — Rate of formation of white dwarfs and Type la super- 
novae, as a function of star formation rate; all rates are normalized 
per unit luminous stellar mass. The left panels are for passive (non- 
star-forming) galaxies (E/SO's) with an age of 8 Gyr (solid line) and 
extending from 5 to 11 Gyr (shaded area). The right panels are 
for objects with star formation. At the bottom of the diagram, 
data points show the observed formation rate per u nit host galaxy 
mass of Type la supernovae from the SNLS survey ISullivan et al.l 
as a function of star formation rate per unit mass. The 
upper shaded area gives the locus of models (with a wide range 
of ages and star formation rate histories) for the formation rate of 
white dwarfs. The solid lines are models with SFR <x , with 
p = 0.95 (rapidly decreasing SFR) on top and p = — 1 (rapidly in- 
creasing SFR) at the bottom. Age runs along each model line, from 
0.1 Gyr at upper right to 10 Gyr at lower left. The star symbols 
are composite stellar populations (old bulge/halo plus star-forming 
disk, with proportions chosen to match observed colours along the 
Hubble sequence (Buzzoni 2005); from left to right the symbols 
represent Sa, Sb, Sc, Sd, and 3 starburst galaxies (p = —0.8) with 
ages 10 Gyr, 3 Gyr, and 1 Gyr. Scaling the formation rate of white 
dwarfs by a factor jj^jg = 0.008 provides an excellent match to the 
observations at all star formation rates (including passive galaxies) . 

IMF. For a lSalpeteij (|1955D mass function dN/dm oc 
(a = —2.35), and a power-law approximation to evolu- 
tionary timescales tevoi ^ lOi°(M/M0)'' yr (6 = -2.5), it 
can be shown that white dwarfs form from a pure instan- 
taneous burst of star formation at a rate that decreases 
with time roughly as ^('i-f'+i)/^' This (somewhat 

counterintuitive) result is due to the fact that the evolu- 
tionary timescales of massive stars are so much shorter 
than for low mass stars, even in the presence of a steeply- 
sloped IMF favouring low mass stars. 

The results of a detailed calculation of WD formation 
rate are shown in the upper pane ls of Fig. 1, using 
more exact evolutionary timesc ales (iBuzzonil jpOna) . and 
a power-law SFR oc t-^ (e.g. IBuzzonil \2mm . (Other 
SFR(t) functions give equivalent results, as will be 
shown.) For passive systems with no star formation (left 
hand panels), we plot the WD formation rate per unit 
galaxy mass (WDR/M) for a star burst with age 7.9 Gyr 
(the approximate age of the Universe at the mean red- 
shift of the SNLS SNe). The shaded region in Fig. 1 
includes systems as young as 5 Gyr (which would have 
formed > 2 Gyr after the big bang for the most distant 
objects in SNLS). For active systems with star formation, 
we plot WDR/M as a function of star formation rate per 
unit galaxy mass (SFR/M). 

Remarkably, the form of this diagram is very insen- 
sitive to the details of star formation history (other 
than age). Old starburst models with increasing SFR(t) 
{p < 0) overlap young models with decreasing SFR 



{p > 0), and all of these models lie within the upper 
shaded region of the figure. A variety of composite stel- 
lar populations, comprising an old passive bulge/halo 
and a star-forming dis k (with the p roportion of bulge 
to disk mass given by Buzzonill2005L to match observed 
galaxy colours) can be seen to fall in the shaded re- 
gion of Fig. 1, as do models with bursts of star for- 
mation. Models with exponentially decreasing SFR(t) 
match the locus of the SFR oc t~P models, as do mod- 
els construct ed using the PEGASE.2 ev olutionary in- 
fall scenario (|Le Borgne fc Rocca-Volmerangc 20021 . Fi- 
nally we note that the range of B-K colours exhibited 
by all of these models spans the full dynamic range of 
colour exhibited by real galaxies. We conclude that the 
locus of the white dwarf formation rate shown in Fig. 1 
is well-determined, independently of uncertainties in the 
evolutionary models. 

SFR/M has units of T~^, and is approximately the in- 
verse of the gas consumption timescale of a galaxy (if 
SFR were held constant). Identifying r ~ (SFR/M)-i 
as a characteristic mean age of a model, one expects 
(WDR/M) ~ ~ (SFR/M)i/2, roughly as observed 

in Fig. 1. The increase of WD formation rate with 
SFR/M is therefore a reflection of the fact that the mean 
age of stellar populations with a large SFR/M (small gas 
consumption timescale) is smaller. 

2.2. Type la Supernova Rates and the Conversion 
Efficiency 

These calculations of specific white dwarf formation 
rate, WDR/M, can now be compared with estimates of 
the SNe la rat e per unit stellar mas s as derived from the 
SNLS survey (jSullivan et al.ll2006f ). Scahng the white 
dwarf formation rate curve by a factor rjobs = 0.008 pro- 
vides an excellent match t o the SNLS observatio ns (and 
also to the observations of lMannucci et ani2005f ). 

However, to convert rjots (the ratio of SN la to WD 
formation rates) to the physically more interesting con- 
version efficiency r] (the fraction''^ of white dwarfs that 
explode as SNe la), it is necessary to correct for the time 
delay, At\YD-SN , between the formation of a WD pri- 
mary and an SN la event. In the case of an SD system, 
this delay is due to the residual time that it takes the sec- 
ondary star (of mass M2 = qMi , where q is mass ratio) 
to evolve. The avail able evidenc e suggests that the distri- 
bution of q is flat (Larson 20( )lD, or even st rongly peaked 
towards equal masses (Pins onneault fc Stanek 2006 1. for 
close binaries. 

We have computed detailed models for SD binaries 
with different distributions of q, using both numerical 
integration and Monte Carlo techniques. A flat distri- 
bution in q results in rates that are < 2x lower than 
for undelayed q = 1 models (see Fig. 2). Most im- 
portant, this factor is nearly constant for all SFR/M. 
Correcting ijobs for the distribution of q values yields 
a white dwarf to SN la conversion efficiency 77 in the 
range 0.01-0.015. (A power-law distribution of mass ra- 
tios N{q) oc q^, corresponding to a distribution of sec- 
ondary masses N{M2) oc Mf , results in constant 77 pro- 

Note that rj includes the binary fraction — e.g. if the fraction of 
binaries were 0.5, then the fraction of WDs in binaries tha t become 
SNe la would be ~0.02. Available evidence IjLadal 120061 ) suggests 
that the fraction of binaries increases with mass. 
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Fig. 2. — Effect of a distribution of secondary masses for the 
SD model. The shaded region is from Fig. 1. The models all have 
fiobs = 0.015 (double that in Fig. 1) and SFR oc t-0-95. The top 
solid line is for binaries with mass ratio q = 1 (as in Fig. 1 - 
no delay between WD and SN). The other two solid lines are for 
populations of binaries with a fiat distribution of q over the range 
[0.5,1] (middle), and [0.2,1] (bottom). The dashed/dotted lines are 
mildly increasing/ decreasing (by a factor of 2) q distributions over 
q =[0.2,1]. It can be seen that the effect of a distribution of mass 
ratios is to decrease SN rates at a given SFR by at most a factor 
of ~2. 

vided that x > —1.5. A distribution of secondary masses 
drawn from tlie IMF gives WDR/M~constant, indepen- 
dent of SFR/M. Th is is discussed further in §3 and in 
iPritchet et all (|2008l Paper 2).) 

A calculation for the DD model is beyond the scope 
of this paper. However we note that a similar result 
(77 ^constant) is obtained if the distribution of semimajor 
axes for DD binaries is independent of mass (Paper 2). 
We also note that 77 « 0.01 is valid for any mass range 
dominated by the SD channel. 

An alternate way of describing the data is through the 
delay time distribution (DTD). From §2, it can be seen 
that the data are consistent with a continuous delay time 
distribution DTD cx (for a burst). Fig. 3 shows that 
the data are bracketed by DTD oc t'"-^ and t'^-'^. 

In summary: The SNIa rate is approximately 1% of 
the stellar death rate, from passive galaxies (SFR=0) to 
starburst galaxies (SFR/M > 10~^ y'^), with httle de- 
pendence on the details of the stellar population mix. 
The dependence of SN la rate on star formation rate in 
both active and passive galaxies is consistent with a sim- 
ple timescale model based on the SD scenario, with a 
uniform WD conversion efhciency rj wl% for all stars. 
The SN la rate - SFR rate dependence is also consistent 
with a continuous DTD oc i-0-5±0-2 jf niore than one 
channel produces SNe la, 77 « 1% still holds for any pro- 
genitor mass range in which the SD channel is dominant. 

Our models are different from those of Greggio (2005) 
since we consider the scaling of the SN la rate with SFR, 
thus minimizing the effects of variations in star formation 
history. In addition our models do not initially include 
variations of 77 with mass (for reasons that will become 
apparent). However, the slope of the DTD for a simple 
aging burst agrees with that of lGreggiol (|2005l . Fig. 3) 
for t <1 Gyr. 

2.3. Uncertainties in the Conversion Efficiency 



What is the uncertainty in the efficiency factor r/? 
Changing the functional form of SFR(t) has already been 
eliminated as a source of error. Changing the IMF to a 
two-piece or three-piece IMF ()Kroupal l2007f ) (with in- 
creasing power-law slope at higher mass) remains con- 
sistent with the locus of white dwarf models in Fig. 1, 
provided 77 is reduced by 25%. Systematic errors in mass 
affect the vertical placement of data points in Fig. 1, 
and these could be a factor of 2 or larger. Systematic 
errors of a factor of 2 in SFR/M affect 77 by about 40%. 
The overall completeness of the S NLS rates, from whic h 
7/ was derived, approaches 100% (| Sullivan et al .11200^ . 
It thus appears that the efficiency factor r] is accurate to 
about a factor of two. Most important, the relative effi- 
ciency factor between galaxies of different SFR/M should 
be unaffected by these systematic errors. 

2.4. Comparison with Other Values of Conversion 
Efficiencies 

Is ?7 = 0.01 reasonable from other considerations? This 
is a difficult que stion to answer. Theoretical models (e.g. 
lYungelsonl[2005l) give efficiencies in the range 0.001-0.01, 
but are subject to a host of uncertainties (e.g. the dis- 
tributions of binary mass ratios and separations, and the 
precise evolutionary scenario for SNe la). For the Milky 
Way, the observed SN la rate (ICappellaro fc Turatto 
200lh and WD formation rate (cf. Phillips! 120021 : iSokei 
2006[) yield 77 ~ 0.1-2%. This number is extremely uncer- 
tain, but it would seem to rule out conversion efficiencies 
larger than a few per cent. 

Other observational determinations of efficiency (from 
the ratio of SNe la t o core collapse SNe, from the A and 
B rate parameters (jScannapieco fc Bildsten|[2005D . and 
from abundances) span a huge range (1-40 %), even ap - 
proaching 100% for intermediate mass stars (jMaod [2007l). 
Part of this discrepancy is due to the different progenitor 
mass range (3-8Mq) used by Maoz; however our deter- 
mination of ?7 = 0.01 will still apply for any mass range 
dominated by the SD channel. To raise the value of 77 to 
0.1 would require 90% SN la incompleteness in the SNLS 
survey, or SFR values that are overestimated by a factor 
of 100. Both of these possibilities are untenable. 

3. IMPLICATIONS FOR SN lA PROGENITORS 

We now turn to relative values of 77 for populations 
with different SFR. The relative WD to SN la conversion 
efficiency 7/ can be measured directly in Fig. 1 for popu- 
lations with different SFR/M, i.e. mean age or mass. In 
particular, it can be seen that the conversion efficiency 
for old stars in passive galaxies is similar to that for the 
most active starburst galaxies. The variation of 77 with 
SFR is smaller than a factor two. 

The SD model results in lower conversion efficiencies rj 
for low mass stars, for a variety of reasons. A passive 10 
Gyr old population is populated with stars with M < 1 
M0, and, because of mass loss on the giant branch, such 
stars produce white dwa rfs of mass 0.5 M© ([Weidemannl 
I2OOOI : iKalirai et al]l2007f l. A binary system initially com- 
prising two 1 M0 stars would therefore not evolve into 
a 1.4 Mq white dwarf unless mass transfer were excep- 
tionally efficient. Below about 2 Mq (turnoff age 1 Gyr), 
a He (rather than C-l-0) white dwarf is the likely out- 
come in a close binary (|Greggio, 2005). and such stars 
do not produce SNe la. The evolutionary time of the 
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Fig. 3. — Data from lSuUivan et al.l II2006I ') is compared with SFR 
oc t~'P models. Solid lines are models from Fig. 1. These models 
correspond to a delay time distribution DTD roughly oc t^^-^ . 
The upper shaded area is bounded by DTDoc t~^-^ (top boundary) 
and t^"'^ (bottom boundary); these models are normalized at high 
SFR. The dashed line represents a p = 0.95 model with the WD-SN 
efficiency at 1 Mq reduced by a factor of 10. The dashed horizontal 
line and shaded area at the lower left shows the predictions of the 
low efficiency model for the passive galaxies. It can be seen that 
models that are normalized correctly for starburst galaxies fail by a 
factor of 7 or more to account for the SN la rate in passive galaxies. 

secondary sets the clock for the SN la explosion, and at 
low mass, some fraction of the secondaries have not yet 
completed their evolut ion. An increasing binary fraction 
with mass (|Ladall200^ lowers the SN la rate at low mass 
relative to high mass. All of these factors will conspire 
to produce lower efficiencies at low SFR/M, relative to 
the most active galaxies. 

The models of iGreggiol (|2005D represent SD efhciency 
at low mass by including a continuous variation of 77 for 
M2 < 2 M0, with ?7iMq/??2Mq < 0.1. The result of such 
an efficiency variation is shown in Fig. 3. Normalizing 
the models to the most active galaxies, the computed 
and observed SN la rates for passive galaxies disagree by 
a factor of 7. If SD progenitor scenarios are less efficient 
at low mass, as is expected, then SD progenitors cannot 
be the only objects producing SNe la; some other class 
of model (perhaps the double degenerate model) must be 
invoked to explain some, and perhaps all, SNe la . A sim- 
ilar co nclu sion has been reached by, f or example. iGreggid 
(|2005( 1 and lYungelson fc Livid (|2000( ). This conclusion is 
also consistent with the observational lack of hydrogen 
in the late-time ejecta of SNe la (Leonard 2007) . 

If secondary star masses for close binaries were drawn 
from the IMF, then close binaries with a high mass pri- 
mary and low mass secondary would be common, and 
the number of SNe la in old progenitors would be en- 
hanced, compensating for the expected lower conversion 



efficiency at M < 2Mq. We reiterate, however, that 
the available evidence supports an IMF for secondaries 
in close binaries that is not drawn from the IMF (e.g. 
lLarsod[200lHPinsonneault fc Stanekll2006[ ). 

The most sub-luminous SNe la (stretch paramete r s < 
0.8) are neglected in the rates (jSullivan et al.l[2006l ). In- 
cluding these objects (mostly in passive galaxies) makes 
the disagreement between observations and predictions 
of the SD scenario even more extreme. A strong type- 
dependent incompleteness (with 90% incompleteness for 
starburst galaxies) would be required to explain the ob- 
servations; there is no evidence for such an extreme effect. 

4. CONCLUSIONS 

The SN la rate is about 1% of the white dwarf for- 
mation rate, for all star formation rates and for a wide 
range of stellar population mixes. The SN la rate is 
also consistent with a continuous delay time distribu- 
tion ^ ^-o.5±o.2^ rpj^g dependence of the SN la rate on 
SFR matches the predictions of a simple timescale model 
based on the single degenerate channel, but only under 
the (unrealistic) assumption that the fraction (1%) of 
white dwarfs exploding as SNe la is independen t of pro- 
genit o r mass. We conclude, as have others (e.g. IGreggid 
120051 : lYungelson fc Livid [2000[ ). that a class of model 
other than the SD model must be invoked to explain 
at least some SNe la. Nevertheless, the conversion ef- 
ficiency of 1% holds for any mass range dominated by 
the SD channel. This conclusion is valid provided that 
secondary stars in close binaries have a distribution of 
mass ratios q that is approximately flat or rising towards 
q — 1- Clearly there is a need for a better understand- 
ing of the the distribution of secondary masses in close 
binaries. 

Finally, it should be noted that the predictions of SN 
la rate in this paper provide a better match to ob- 
servations tlia2_does_tlie_em£i™ M + B ■ SFR 
formu la (jScannapieco fc Bildsted l2005t ISullivan et al.l 
I2OO6I ). The implications of this will be discussed in Paper 
2. 
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